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Abstract 

Point-to-point reflection holding for harmonic functions subject to 
(N| the Dirichlet or Neumann conditions on an analytic curve in the plane 

almost always fails for solutions to more general elliptic equations. 
We develop a non-local, point-to-compact set, formula for reflecting 
a solution of an analytic elliptic partial differential equation across a 
real- analytic curve on which it satisfies the Dirichlet conditions. We 
also discuss the special cases when the formula reduces to the point- 
T"! to-point forms. 

> 

X 

1 Introduction 



> 
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ON 



Schwarz symmetry principle is one of the celebrated tools in analysis and 
mathematical physics that has been attracting attention of many mathemati- 
cians pQ- [2], [H] [20], [22] [2E]- From the point of view of applications it 
is important to have an explicit reflection formula for a specific problem ([7], 
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|lUj . [22]). One of the open questions is the following: for what partial differ- 
ential equations, boundary conditions and spatial dimensions such a formula 
exists and what is the structure of this formula, in other words, whether it is 
a point to point formula (see, for example [8J) or it has a more complicated 
structure, for example, a point to a finite set [20J or a point to a continuous 
set (see, for example, [2] and references therein). 

In this paper, we derive a reflection formula for solutions of elliptic equa- 
tions in IR 2 with respect to a non-singular real-analytic curve and study the 
obtained formula. We call this formula non-local, since unlike the classical 
point to point reflection (see the Theorem |1.1| below) this is a point to com- 
pact set reflection, generalizing the following celebrated Schwarz reflection 
principle for harmonic functions. 

Theorem 1.1 fJT7|/ Chapter 9, p. 51; [28] Chapter 1, p. 4). Let T = {(x,y) : 
f(x,y) = 0} C IR 2 be a non-singular real-analytic curve and P' G T. Then, 
there exists a neighborhood U of P' and an anti-conformal mapping R : U ^ 
U which is identity on T, permutes the components Ui, U2 ofU\T and relative 
to which any harmonic function u(x,y) defined near V and vanishing on F 
is odd; i.e., 

u(x ,y ) = -u(R(x ,y )) (1.1) 

for any point (xq, yo) sufficiently close to T. Note that if the point (xq, yo) G 
U\, then the "reflected" point R(xo,yo) G Ui. 

Here the mapping R can be described by considering a complex domain 
Uc in the space C 2 , such that [/cflR 2 = U, to which the function /, defining 
the curve T, is continued analytically. After the transformation of the vari- 
ables, z = x + iy, ( = x — iy, the equation of the complexified curve T c can 
be rewritten in the form 



If grad f(x,y) 7^ on T, (1.2) can be solved with respect to z or (; the 

corresponding solutions we denote as ( = S(z) and z = S(Q. The function 
S(z) is called the Schwarz function of the curve T [B], Chapter 5, p. 21. The 
mapping R is given by 



R(x,y)=R(z) = S(z). (1.3) 
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Formula (1.1) has been generalized to cover several other situations in- 
cluding the Helmholtz equation and wave equation, and the polyharmonic 
functions (see, for example, [25], [I], [20] and references therein). The pur- 
pose of this paper is to obtain an explicit reflection formula for solutions to 
the elliptic equation 

du du 

Lu = A X;V u + a(x, y)g^ + K x > y ^Qy~ + c ( x ' y ^ u = ( 1>4 ) 

with respect to a real analytic curve in IR 2 , where the solution vanishes, and 
to investigate the properties of the mapping induced by this formula. Here 
a(x,y), b(x,y) and c(x,y) are real-analytic functions in the domain U C M 2 . 

In what follows, a formula, expressing the value of a function u(x, y) at 
an arbitrary point (xq, yo) G U\ in terms of its values at points in f/ 2 , is called 
a reflection formula. It is more often an integro-differential operator than a 



point-to-point reflection (1.1), which seems to be quite rare for solutions of 
partial differential equations. In particular, for solutions of the Helmholtz 
equation (A x>y + \ 2 )u(x, y) = vanishing on a curve T, point-to-point reflec- 
tion holds only when T is a line, while for harmonic functions in M. 3 it holds 
only when T is either a plane or a sphere [3], [18j. The paper by P. Ebenfelt 
and D. Khavinson [8] is devoted to further study of point-to-point reflection 
for harmonic functions. There it was shown that point-to-point reflection 
in the sense of the Schwarz reflection principle for n > 2 is very rare in R n 
when n is even, and that it never holds when n is odd, unless T is a sphere 
or a hyperplane. Reflection properties of solutions of the Helmholtz equation 
have also been considered in [9], [23], |25j . Two later papers are devoted to 
derivation of non-local formulas for Helmholtz equation subject to Dirich- 
let and Neumann conditions respectively. Recently a reflection formula for 
harmonic functions subject to the Robin condition, ad n u + (3u = 0, on a 
real-analytic curve was derived in [2], and it was shown that the obtained 
(non-local) formula reduces to well-known point-to-point reflection laws cor- 
responding to the Dirichlet and Neumann boundary conditions when one of 
the coefficients, a or /3, vanishes. 

The structure of the paper is as follows: in Section [2] we describe some 
preliminaries; in Section [3] we formulate the main theorem, which is proven 
in Section |4} Conclusions and the special cases, when the point to point 
reflections hold, are discussed in Section [5j 
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2 Preliminaries 



We are starting this section by recalling a classical B. Riemann's result for 
hyperbolic equations (see [13] Chapter 2, p. 65 or [T2] Chapter 4, p. 127 for 
detailed explanations, here we follow a short version [T7| Chapter 9, p. 55). 
Consider a hyperbolic differential equation with entire coefficients 



d^u du du 

HU ^ dxdy + a(yX ' V ^dx + V " > dy + V ^ U = °' 



its adjoint equation is 



H*u 



d 2 u d(au) d(bu) 



+ cu = 0. 



(2.1) 



(2.2) 



dxdy dx dy 

The Riemann function 91h(x, y; xo, yo) of operator H is defined as the solution 
to the Goursat problem: 



H*9K H = near (x ,y ), 

f v 

9i H (x , y, Xo, yo) = expj / a(x , r)dr 

yo 

M H {x,yo;x ,y ) = expjjf 6(i,y )dt}. 



(2.3) 



Note that %\h is an entire function of all four variables, moreover 9t H (x, y; x , y ) 
y^H*(xo,yo] x,y), 9tH{xo,yo',xo,yo) = 1 and the following Riemann's lemma 
holds. 

Lemma 2.1 Let T := {(x,y)\ y = s(x)} be a non- characteristic with respect 
to H real- analytic curve that divides a domain U C 1R 2 in two connected 
components U\ and U2, and u(x,y) be a solution of (2.1) near T. For all 
points P(x ,y ) G U sufficiently close to Y we have 



N 



u(P) = ^u(M)K H (M) + ^u(N)m H (N) - J (IXdy - 2Jcfe), 



(2.4) 



where M 



1 (yo),yo), N = (x ,s(x )) and 



^ 1 du ld<R H 
2 ay 2 ay 

V3 = bm H + -m H —--—^u. 

2 dx 2 dx 
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If in addition solution to the equation Hu = vanishes on T, formula 
u(P) 



(2.4) reduces to 



N 



- 2 j^n{ Yx dx--dy). 

M 



(2.5) 



Remark 2.2 For i/ie wave equation, a = b = c = m (2.1), the Riemann 



function equals 1 identically. Consider a point P(x , y ) G U\ and a solution 
of the wave equation vanishing on V. Let's allow the path of integration 



MN in (2.5) to degenerate to a pair of segments (in U2) of a vertical and 
horizontal characteristics through points M and N, which intersect at a point 
Q(s^ 1 (y ), s(xq)) G Ui. Then formula (2.5) becomes 



u{P) = -u{Q). 

Since the points P and Q are located on the opposite sides of the curve T, 
the later formula states a point to point reflection law for the wave equation. 

Remark 2.3 If for a solution of the wave equation vanishing on T we allow 



the path of integration MN in (2.5) to degenerate to a polygonal line consist- 
ing of vertical and horizontal segments with vertices M = Qi, Q2, Q3, . . . , Q n = 
N such that Q2k+\ G V and Q2k G U2 \ T, k = 1, . . . , {n — 1)/2, then a version 
of a point to finite set reflection will be obtained [20] 



u{P) 



(n-l)/2 



2k I 



k=l 



(see JW^ for other examples of point to finite set formulas). 



If we consider the elliptic equation (1.4) in the complex domain Uc C C 



then the equation and its adjoint in characteristic variables (z, Q become 



similar to the hyperbolic equation (2.1) and its adjoint (2.2) 

du 



L — + A^ U + P 

dzdC, dz d( 



Cu = 0, 



, : _ d 2 u d{Au) d{Bu) 
L r u = - - ; — h Cu = U, 



dzd( 



dz 



(2.6) 
(2.7) 
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where the coefficients in (1.4) are replaced with 



I i 

A{z, C) = ^ [a(x, y) + ib(x, y)] , B(z, C) = - [a(x, y) - z&(ac, y)j , 

C(z,C) = ^c(x,y). 

Analogously, the Riemann function of L is defined as the solution to the 
Goursat problem in C 2 : 

< Li*R=£- c tR-£ z (AM)-£ c (Bm) + cm = o, 



<R\ Z=ZQ =exp{/A(2!o,r)tir} J ^ 
$H, c=Co =exp{/5(t,Co)dt}. 



By a fundamental solution of operator L we understand a solution of the 
equation L*G(xo, yo, ^ ? y) = 8(%o, Vo), where L* is the adjoint to L differential 
operator. Thus, function G written in the characteristic variables z = x + iy 
and C = x — iy is a solution to the equation 

d 2 u dAu dBu _ ^ . . 

Llu = lMK-^-lK +Gn = 5M) - (2 ' 9) 

The following formula (see [13] Chapter 3, p. 72) shows that the Riemann 
function is a factor of the logarithm in an expression for the fundamental 
solution of the operator L c 

G(z,C;*b,Co) = -^(^C;^o,Co)ln[^-^o)(C-Co)]+^,C^o,Co), (2.10) 

where t? d z o> Co) is an entire function. 

Note, that the fundamental solution exists (see [TB] Chapter 3, p. 50) and 
is uniquely determined up to the kernel of operator L. 

There are different representations of the fundamental solution, for ex- 
ample, [5]; [15] Chapter 3, p. 76; [IB] . However, for what follows we need 
a special representation as a sum of two functions, each of those has the 
logarithmic singularity on a single characteristic in C 2 . This representation 
is given by the following theorem. 
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Theorem 2.4 There exist a fundamental solution of L that can be rep- 
resented in the form 

1 ~* (2.11) 

(2-12) 



g — ———{G\ + cy, 

An 



G o = ^2 a i( x o^yo : > x ,y)fk( , >Pj), j = 1,2 



k=0 



A(0 



(_l)-fc-i(_jfc _ i)!^ ; fc<-l, 



^(lne-Cfe), fc = 0,l,..., 
C = 0, C k = J2\i k = l,2,... 



(2-13) 



^ = (x-x )+i(y-y ) = z-zo, ip 2 = (x-x )-i(y-y ) = C~ Co- (2-14) 

i/ere i/ie coefficients a 3 k are uniquely determined by recursive transport equa- 
tions 



£a4 = 0, ^ a k+i = ~~ L*cO? 



k ' 



d 



dz ld( 
subject to the initial conditions 



A 



d_ 
.dz 



B 



(2.15) 



a 



°lc=c 



exp|/5(t,Co)d*}, "o 



0, k = 1,2, 



-0 

C 



expj / A(z ,T)dr\, 



°lc=c 

^|_ =0, fc = l,2, 



(2.16) 



Note that (2.15) and (2.16), in particular, imply 

C C 
«o = exp (^J A(z, r)dr+ J B(t, Co)dij , a 2 Q = exp ( / A(z , r)rfr+ ^ £(t, C)dfc) • 



Co 



20 



Co 



20 



(2.17) 

Taking into account (2.10 ), one can interpret a J k as coefficients in the following 
series representations for the Riemann function (2.8) [26J: 

00 ( z — z )k °° ({ (* \k 

W-(zo,(o,z,() = ^2 oil(z , Co, z,() = ^ a 2 k (z , (o, z,() Tjp— ■ 

fc=0 ' A;=0 

(2.18) 
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Remark 2.5 For the Laplace equation, a = b = c = 0, and, therefore, 
A = B = C = 0. Thus, al = a% = 1 and a] = a| = 0, j > 1, and 

Gf = ln(*-*b) and G\ = ln(( - Co) (2.19) 
respectively, which leads to a standard fundamental solution: 

G=~ ln[(* - «b)(C - Co)] = ln[(x - x ) 2 + (y - y ) 2 )]. (2.20) 

Remark 2.6 In £/ie case of the Helmholtz equation a = b = ana 1 c = A 2 ; 
t/jat is ; ■^q^u h + = 0. Here \ is a real number, functions G% and G 2 

reduce to the form used in WM 

Summing up Gf and G^ and multiplying by — ^ one obtains well-known 
fundamental solution of the Helmholtz equation: 

-^(Gf +Gf ) = C + ^ V2 Jo(A^ - * )(C ~ Co))-\N (\y/(z - z )(C - Co)), 

(2.22) 

where c is the Euler constant, and Jo and iVo are the Bessel and the Neumann 
functions of zero order respectively. 



Ck 



(2-21) 



3 The main result 

Consider a solution of homogeneous linear elliptic differential equation, writ- 
ten in its canonical form [12] Chapter 5, p. 136 (with the Laplace operator, 
Ajbj,, in the principal part), in a domain [/cl 2 vanishing on an algebraic 
curve T, 

Lu = A x<y u + a d f x + b d f y + cu = near T, 

t/)| r = 0; a, 6, c are real-analytic functions of x, y. 
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Theorem 3.1 Under the above assumptions, the following reflection formula 
holds in U: 



u(P) = -c (P,I>(Q) + 



2i J y- ox ox * oy oy * 



(3.2) 



where P = (xo,yo) and Q = R(P) (see (1.3)), and the integral is computed 
along any curve joining V with Q. Here 



c (P,r) 



cxp 



+ exp 



B(t, S{z ))dt + J A(z ,r)dT 

z o Co 
S(z ) S(Co) 

A(S(( ), r)dr+ [ B(t,(o)dt 



(3.3) 



Co 



~o 



where 



A ( z > = 4 [ a ( x , y) + ib ( x , y)] , B ( z > = \ [ a ( x > y) - ib ( x , v)] ■ 



C{z,C) = jc(x,y), V = V(x ,y ,x,y) = V 1 (x ,y Q ,x,y) - V 2 {x ,y ,x,y). 
Functions Vj are solutions of the Cauchy-Goursat problems: 

J = 1,2, 
J = 1,2, 



V x = expl / A(S(C),r)dr + / B(t,()dt\, on the char. l x = {S(() = z } 





= 0, 








C 


, *i= 


expj J 




Co 




C 


v 2 = 


expj J 




Co 



V 2 = exp< / A(z,r)dT+ / B(t, S{z))dt >, on the char. l 2 = {S(z) = Co} 



20 



where is the adjoint operator to Lc and 9i(zo, Co, z, C) the Riemann 
function of L. 
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4 Proof of the Theorem 



3.1 



4.1 Sketch of the proof 

We begin with Green's formula expressing a solution of the equation Lu = 
at a point P via its values on a contour 7 C U\ surrounding the point P |TT] : 



u{P)= / u)[u,G], (4.1) 



where 

u[u, G] = {u^- - G^- - auG] dy - {u^- - G^- - buG] dx. (4.2) 
ox ox oy oy 

Here G = G(x,y,xo,yo) is an arbitrary fundamental solution of L, that is, 
a solution to the equation L*G(xo,yo,x,y) = 5(xo,yo). It is well-known 
that G is a real-analytic function in R 2 except at the point P(xo,yo)- Its 
continuation to the complex space C 2 has logarithmic singularities on the 
complex characteristics passing through this point, i.e., on K P := {(x — x ) 2 + 
(y ~ Ho) 2 — 0}. Our proof is based on the idea suggested by Garabedian [TT] 
to deform contour 7 across the curve T from the domain U\ to the domain U-i- 
To be able to realize this deformation, first, we use a special representation 
for a fundamental solution, that is, a sum of two functions, each of those has 
a singularity on a single characteristic only. This representation is given by 



the Theorem 2.4 above. Next, we replace the fundamental solution G with 
a so-called reflected fundamental solution. After proving the existence and 
uniqueness of the reflected fundamental solution, we describe the deformation 
of 7 and obtain the desired reflected formula. Finally, we simplify the formula 
and discuss the cases for which it reduces to the simplest point to point form. 



4.2 The reflected fundamental solution 

This section is devoted to the construction of the reflected fundamental so- 
lution G, which plays the key role by enabling us to deform the contour 7 
across the boundary. G depends on the operator L and the curve r ^\ As 
it will be shown in the next two sections, the reflected fundamental solution 

1 G depends on the boundary condition as well, but the later is beyond the scope of 
this paper, see [2j, [25] and [27] for some relevant results. 
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determines whether the corresponding reflection formula can be reduced to 
the point to point form. 

Function G is a solution of the equation L* C G = subject to the boundary 
condition G = G on Tq and has singularities only on the "reflected" charac- 
teristic lines l\ and 1% (see Fig. 1) intersecting the real space at the reflected 
point Q = R(P) in the domain U 2 and intersecting Tq at Kp fl Tq- 

We seek the reflected fundamental solution in the form 

G(zq, Co, z, C) = - — (Gi(zo, Co, z, C) + G 2 (z , Co, z, ()), (4.3) 

where the functions Gj, j = 1, 2 are defined as the solutions to the following 
Cauchy-Goursat problems with prescribed singularities, 



L* c Gj — 0, 



1,2, 



G; 



Gq, 



(4.4) 



Gj has singularities only on the char. lj =: {ipj{z, C) = 0}, 



where ipi = S(()—z and ip 2 = S(z) — Co are solutions of the Hamilton- Jacobi 
equation 

dt/jj d'ipj 



dz d( 



0. 



(4.5) 



First, we construct the solutions to the problems (4.4) as some formal 
expansions. Then we justify their convergence. 
We seek this expansion in the form 



Gi = /3l(z 0} (o, z, ()f k(ipi] 



k=0 

oo 



G 2 = ^2 Pkizo, Co, z, OfkM, 



(4.6) 
(4.7) 



k=0 



where /&(£) is defined by (|2.13[). 



Substituting (4.6) and (4.7) into (4.4), we obtain the following recursion 



for the coefficients f3\ and 

Wo 

dz 

aft 



BPl)S' (C) = 0, 
Affi)s'(z) = 0, 



d_ 

dz 



Pl +1 -BPl +1 )S'(() = -L* c /3l k>0 
0l +l -A(3l +1 )s'(z) = -L* c /3l 



k > 
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subject to the following initial conditions 



PL 



a 



a 



k = 0,1,2,... . 



(4.9) 



Note that both S'(z) and S'(C) do not vanish on T c , [5] Chapter 7, p . 42. 
Indeed, functions S(z) and S(C) are inverse of each other (see (1.2)), so 
S(S(z)) = z. Differentiating the later equation and taking into account that 
S(z) = C on r c , we obtain S'(C) • S'{z) = 1. In R 2 , S'(£) =jS'(z), therefore, 
= |S"(C)I = 1 on T. Thus, both functions S'(z) and S"(C) are nonzero 
throughout some neighborhood of T as continuous functions. 

Thus, functions j3\ and /3| are uniquely determined near Tq, specifically 

C z 5(C) 

ft = exp([ A(S(C),r)dr+ [ B{t,()dt + / [B(t, Co) - B(t, ()}dt 



Co 



20 



2(1 



C S(z) 

/?o=exp^y A(2,r)rfr+ ^ [A(z , r) - r)]dr + ^ B(t, S(z))dtj . 



Co 



Co 



20 



(4.10) 

Hence, the formal expansions for the functions G\, G<i satisfying conditions 
(4.4) are constructed. 



Lemma 4.1 The series (4.6) and (4.7) converge near I\ 
Proof of the Lemma 14.11 



Let us prove the convergence of the series (4.7) by considering an auxiliary 



family of problems depending on parameter £ : 
L* c Vt(z o ,Co,z,C,Z) = 0, 

Vz(zo,(o,S((o-£),C,£) = 0. 



(4.11) 



Here $ is a given analytic function, that has Taylor expansion 

(S(z)- Co + fc+1 

k=0 



$(z ,(o,z,S(z),£) = ^Tal(z ,(o,z,S(z))- 



(k + l)\ 



where coefficients aq.(zo, Co, z, C) are the same as in (2.18) [26] . 
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Taylor expansion of the solution to the problem (4.11) (if it exists) has 
the form 



Vz(zq, Co, z, C, = ^2 £°» z ' 



fc=0 



(S(z) - Co + Q 

{k + iy. 



k+1 



(4.12) 



where the coefficients /3| are the same as the coefficients in series (4.7). Con- 
vergence of the later, therefore, followed from convergence (4.12). To show 



existence and uniqueness of the solution to the problem (4.11) in the class of 
analytic functions we use the substitution 



Vt(zo,Co,z,C,Z)= J dr J n(z ,(oit,T,Odt+$(zoXo,z,S(z),0 (4.13) 



with unknown density /i, which reduces the problem (4.11) to the Volterra 
integral equation 



K z o,(o,z,(,0 + A (z,()S'(z) J n(z Q ,(o,t,S{z),£)dt 

s(&-e) 

c 

-A( z >0 J M^o,Co, z,r,^)dr - B(z,() J n(za,(o,t,(,Qdt (414) 



S(z) 

c 



S(Co-0 



F{z,() J dr J ii(z ,(a,t,T,£,)dT = ^(z ,(o,z, (,£,), 



S(z) S(Co-fl 



where 



F(z, C) = C) + ^B(z, C) - C7(z, C), (4.15) 

(*o, Co, 2, C, = ^(2, C)$(*o, Co, *, + C) JUfo, Co, 2, 0- (4-16) 



The existence and uniqueness of the analytic solution of equation (4.14) can 



be proven by iteration technique described in Chapter 1, p. 11. Thus, 



there exists unique solution of (4.11 ), which has unique Taylor expansion with 
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respect to variable £ at the point £ = —(S(z) — Co), this expansion coincide 
with the expansion (4.12). Thus, series fl4.12[ ) converges in the neighborhood 
of T, and so does (4.7). 

Analogously, considering the following auxiliary problem depending on 
parameter r\: 



~ <x> ~- 

V v (z ,(a,S((),(,v) = E a l( z o, Co, S '(C),C) 

fc=0 

V v (zq, Co, 2, 5(^0 - ??)), ^) = 0, 



(S(C)-*o+»?) fc+1 

(k+iy. 



(4.17) 



whose solution has Taylor expansion V v = EaUo 0k ^%fc+i)^ + • 
show convergence of (4.6). That finishes the proof. 



one can 



4.3 The reflected fundamental solution as a multiple- 
valued function 

As it was conjectured in [I]: "Perhaps Looking-glass milk isn't good to 
drink" . In this section we show that the reflected fundamental solution (the 
looking-glass fundamental solution), except for some special cases, does not 
inherit all of the properties of a "true" fundamental solution, in particular, 



the representation (2.10) with the Riemann function as a factor of the log- 
arithm does not hold. Moreover, as w e are abo ut to show, the factors of 
the logarithms in G\ and G<i (see (4.6) and (4.7)) are not the same, which 



ution a multiple- valued function even in 



makes the reflected fundamental so 
M 2 . The later explains (see Section 4.4) why point to point reflection almost 
always fails. 

Indeed, consider a point moving along a continuous curve 7 surrounding 
either the branch line z = Zq of G\ or the branch line C = Co of Gi (2.11 ). As 



the point makes a complete cycle around the line, it passes to the next sheet 
of Riemann surface, while going around the cyclic path surrounding both 
characteristics at once, it remains on the same sheet of the Riemann surface. 
For the reflected fundamental solution G, the point passes to the next sheet 
even if the curve 7 lays in M 2 and surrounds both intersecting branch lines, 
S(() = zq and S(z) = Co- 

To show this, let us compute the increment of the function G when a curve 



14 



7 C M 2 is a circle of a small radius p centered at the point Q(S(zq), S(z )y. 



oo 
J=0 



i (5(C) - *b) 



(S(z) 



(4.18) 



Taking into account that ( = z in IR 2 , let us set z = S(z ) + pe 1 ^ and 
C = S(zq) +pe~ 11 ^, and expand the Shwarz function and its inverse into Taylor 
series at the point Q: S(z) = zq + Cipe 1 ^ + o(p), S(z) = zq + Cxpe" 1 ^ + o(p). 
Without loss of generality assume that the coefficients A, B and C in 



(2.6) are constants (otherwise we should use their Taylor expansions in this 



analyses), then 0$ = (3q (see (4.10)), and 

ft = {AB - C )e^-^B(z- Z0 ))^ z _ 5(C))/5'(C) + C - Co), (4.19) 
ft = (AB - C) e (A(C-Co)+B(^-^o)) ^ _ S(z))/S'{z) +z-Zq). (4.20) 



Thus, the increment (4.18) becomes 
P 



2i 



(AB - C) [dSizo) - C lZo + S(z ) - z ]C ie 



-i<f> 



■[CtSizo) - C lZo + S(z ) - z }C ie ^ ) + o(p). 



(4.21) 



Formula (4.21 ) shows that the increment can be equal zero only in two cases: 
either when (i) AB — C = or (ii) expressions in the brackets equal zero. 
The later happens if boundary T is a segment of a straight line, while (i), for 
example, holds if operator L is the Laplacian. 

Having the detailed description of the reflected fundamental solution we 
are ready to derive the reflection formula by explaining how the contour 7 in 



(4.1) can be deformed from one side of the reflecting surface Tc to the other. 



4.4 Deformation of the contour 



Formula (4.1) involves integration over a contour 7 C U\ surrounding both 
characteristics on which functions G\ and G2 have singularities (lines l\ and 
I2 Fig. 1). To express value u(P) in terms of values of u(x,y) in U2, that is, 
to construct a reflection formula, it is sufficient to deform the contour 7 from 
the domain U\ to the domain {7 2 - Note that since the integrand in (4.1) is a 
closed form, du = 0, the value of the integral will not change while we are 
deforming the contour 7 homotopically. 
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Figure 1: Contour deformation 



First, the contour is deformed to the complexified curve Tq. Taking into 
account that the characteristics of G passing through the point P intersect 
Tq at two different points in C 2 , assume that the point P lies so close to the 
curve r that there exists a connected, univalently projected onto a plane, 
domain Q C that contains both points of intersections 



We start the deformation with stretching the contour 7 (see (4.1)) in the 



real plane until its small arc reaches the curve T (it becomes a mirror image 
of 7 in Fig. 1). Then we substitute a sum of G\ and G 2 for G in (4.1) and 
split the integral: 



u(P) 



u)[u, Gi] + / cj[u, G 2 ] 



(4.22) 



Note that contour 7 is not closed on Riemann surfaces of each G\ and G 2 



(see Section 4.3). As a point of disconnection (one of two endpoints) let us 
choose a point K e 7 fl V. Then in the first integral in (4.22) we "lift" the 
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contour 7 to 7' (solid line above the plane in Fig. 1) such that we do not 
move some points of 7 D T in the neighborhood of the point K. Then we do 
the symmetric (with respect to plane M 2 ) deformation in the second integral 



in (4.22) 



Taking into account that U\ T = 0, differential form u (4.2) on Tq becomes 



du 



du 



oj'[u, Gj] = Gj — dx — Gj — dy, j = 1,2. 



3 dy ^ ^ J dx 



(4.23) 



(4.4) 



Now we can replace Gj with Gj (see formula (4.3)). Indeed, according to 

u/[u,Gi] = J u/[u,Gj\. (4.24) 

In order to deform contour 7' from Tc to the domain U2, it is necessary to ap- 
ply the "mirror" deformation procedure. Note that during this deformation 
the point K is fixed and contour surrounds one of the "reflected" character- 
istic lines l\ or l 2 (see Fig. 1) intersecting the real space at the reflected point 
Q = R(P) in the domain U 2 and intersecting r c at K P fl T c . 
Finally, we have 



u{P) = / w[u,Gi] + / u[u,G 2 }. 



This formula can be rewritten as a single integral 



u(P) = j u{u,G}, 



(4.25) 



(4.26) 



but as it was discussed in Section 143] the contour 7, generally, is not closed 
on Riemann surface of G, so in most of the cases we do not expect to be 



able to move the point K (see Fig. 2) from the curve Y. Formula (4.26) is a 



version of a desired reflection formula. In the next section we simplify it and 
show that it holds in the large. 

4.5 The reflection formula in the large 



Formula (4.26) in (z,Q variables has the form 

u(P)= [ u>[u,G\, 



(4.27) 
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where 



u[u,G]=i({u^-G^-2AuG}d(-{u^-G^-2BuG}dz^. (4.28) 

Here G is the reflected fundamental solution and contour 7 C U2 surrounds 
the point Q (see Fig. 2). Recall that G is a sum of two series (with certain 




Figure 2: Contour transformation 

radii of convergence). Now we are going to show that the formula holds in 
the large. 

Let us rewrite functions Gi in the form: 

G^Vtln^ + Vt, 1 = 1,2, (4.29) 

where _ _ 

V, = ±?'M V-t0C„ (4.30) 

and 

^ = S(C) - z , ^2 = 3{z) - Co- (4.31) 
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Substituting Q4.29[ ) and ( |4.30[ ) into ( |4.27[ ) and letting the radius of the arc NM 
to zero (see Fig. 2) result in vanishing integrals of the terms involving prod- 
ucts of function V\ and derivatives of function u as integrals of holomorphic 
functions over a closed contour. Combining terms in (4.27) involving deriva- 

ithmic 

(4.32) 



tives of logarithms and separating them from the terms involving logarithmic 
functions yields 

u(P) = Q + I, 



where Q and I in characteristic variables have the form 

E / (uV^iln^dC-uV^ln^dz). 



i 

An 



(4.33) 



I = -- 



An 



oj[u, G\] In ipi 



1 

An 



oj[u, G 2 ] ln-^2 



(4.34) 



Substituting series (4.30) for Vi into (4.33) and computing the residues at 
the point Q where the integrand has the simple pole, we have 

S(Co) S(z ) 

n,(0) fexnf 



a Co) oiz ; 

^t(Q)(exp( J B{t,S{z ))dt+ J A{z ,T)dr) 



Co 



S(z ) S(Co) 

+ exp( / A{S{( ), r)dr+ [ B(t,(o)dt 



(4.35) 



Co 



2() 



which holds in the large. 

Using properties of the logarithmic function and replacing the contour 7 
with a segment EQ, second integral can be rewritten as 



i = 2™(-i-)( / r,j 



Q 



uj[u, V 2 ] 



1 

27 



(4.36) 



E 



where V = V\ — V 2 - Note, that the logarithms in (4.34) have complex con- 
jugated arguments (4.31) in M 2 , however they cancel each other only if the 
factors V\ and V 2 are equal, which generally is not the case. 
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Even though the later formula involves series V\ and Vg, it is also holds 
in the large, since these expansions can be interpreted as solutions of the 
following Cauchy problems 



L* c Vj = 0, 



J = 1,2, 



9\ 



\r c : 
C 



Vj = exp< f A(9j,r)dr + f B(t,r]j)dt }> on the characteristic lj =: {i/ij(z, () 
Co 20 



where 6*i = S((), 2 = z i Vi — C an d V2 = S(z). Problem (4.37) by a 
substitution with unknown density /i, for example, for j = 2 



V 2 {z , Co, z, C) 



C 

dr 



t*(zo,(o>t,T)dt 



S(z) 5(Co) 



(4.38) 



/ A(z,T)dr- / AO,t)o!t 

+ ^(zo,Co,^,^))e c o c ° 



can be reduced to the Volterra integral equation, whose solution as a function 
of four complex variables exists and unique in some cylindrical domain near 
T (see 



Chapter 1, p. 11). Thus, the solutions of (4.37) exist in C 4 as 
multiple-valued analytic functions, whose singularities coincide with those of 
S(z) and S((). 

Combining (4.32), (4.35) and (4.36) we arrive at the formula (3.2), which 
proves the theorem. 



5 Conclusions and remarks 

5.1 Equations with constant coefficients 

We have obtained a reflection formula for elliptic equations with analytic 
coefficients subject to homogeneous Dirichlet conditions on a real analytic 
curve. This is a point to compact set reflection, which in some cases can be 
essentially simplified. 



Consider the case when the coefficients a, b and c in the equation (1.4) 
are constants, 

A XjV u + a— + b— + cu = 0, (5.1) 
ox ay 
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and, therefore, A, B and C are constants as well. In this case solutions oP k 
to the problems (2.15) - (2.16) can be written explicitly, and the Riemann 
function (2.18) has the form 

({z - Zo)(C - Co)(AB - C)) 1 



k=0 



(k\) 



exp(A(C-Co)+S(z-2o)). 

(5.2) 



Our main conclusion confirms the fact that the point to point reflection 
is quite rare. 



Theorem 5.1 For non-trivial solutions of elliptic equation (5.1) with con- 
stant coefficients vanishing on a real- analytic curve T, there is no point to 
point reflection unless one of the following conditions hold: 

(i) T is a line, 

(ii) a 2 + b 2 -4c = 0. 

Proof: The proof immediately follows from the fact that the integral term 
I ^ in fl43~6~| ). Indeed, formula p~21~j ) imply that V ^ 0. Thus, for I to 
be zero, function u and its first derivative must vanish on a path joining the 
curve T with the reflected point, which contradicts the assumption that u is 
not equal zero identically. 

Theorem 5.2 Let V := {ax + (3y + 5 = 0} be a line. Then for any solution 
of the equation Au + au x + hu y + cu = with constant coefficients vanishing 
on T the following point to point reflection formula holds in M. 2 : 

(ax + (3y Q + 5) (aa + b(3) ' 



u(P) 



exp 



u(Q). 



(5.3) 



a 2 + (3 2 

Proof: Under the assumptions of the theorem, the Shwarz function is 
S(z) = mz + q, where 

P 2 -a 2 + i2aP -2a5 + i2fi8 



m 



a 2 + (3 2 



Q 



a 2 + (3 2 



Functions V\ and V2 are equal (see (4.37)), 

((mz + q-( )(rh( + q- z )(AB - C)f 



V 1 = V 2 = J2 

k=0 

and therefore, V 



,(A(C-Co)+B(*-«,)) 



V 1 -V 2 



0, and the integral 1 = (see (4.36)). Formula 
fl4~35) can be simplified, and Q = -u(Q) e Mmz +g-c )+B(fhC+g-z ) \ The 

in variables (x,y) gives (5.3). 
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Corollary 5.3 Let V be a line with equation y = 0, then for any solution of 
(5.1) vanishing on F the following reflection formula holds 

u(x , y ) = -e~ byo u(x , -y ). (5.4) 

Corollary 5.4 Let T be a line with equation x = 0, then for any solution of 
(5.1) vanishing on T reflection formula has the form 

u(x , y ) = -e- axo u(-x , y ). (5.5) 



Corollary 5.5 If a = b = formula (5.3) recovers known point to point 
reflection for solutions of the Helmholtz equation vanishing on a line 

u{P) = -u{Q). 

Remark 5.6 Note, that in the case of the Helmholtz equation, a = b = 
and c = X 2 , when T is a real-analytic curve, formula (4.35) can be simplified, 
and Q = —u(Q), but I ^ in (4.36) unless Y is a line f77[ / Chapter 9, p. 59; 

W; m 



Theorem 5.7 Let T be a real-analytic curve. Then for any solution of the 
equation Au + au x + bu y + (a 2 + b 2 )/4u = vanishing on T the following 
point to point reflection formula holds in M. 2 : 

U (P) = - e ^(*o)-Co)+B(S(Co)- Z o) w (Q). ( 56) 

Proof: In characteristic variables condition c = (a 2 + 6 2 )/4 is equivalent to 
AB — C = 0. Then the Riemann function (5.2) has the simplest form 

^ ,Co,^C) = e A(c - Co)+B(z - Zo) , (5.7) 

and V\ = Vi = 91 for any analytic curve T. Thus, the reflection formula has 
the point to point form (5.6). 

Remark 5.8 Equation Au+au x + bu y + (a 2 + b 2 )/4u = can be transformed 
into the Laplace equation using the substitution u(x,y) = v(x,y)e^^ ax+by ^ 2 , 
where v is a harmonic function, and, therefore, v enjoys the celebrated Schwarz 
symmetry principle (1.1). 



Example 5.9 Formula (5.6) for the unit circle centered at the origin can be 
rewritten in (x, y) variables as follows 

i \ f2{ax + by ){l-x 2 -y 2 )\ x y 

u{x ,y ) = -exp ^— kt — ? ). 5.8 

V x 2 + y 2 J x 2 + y 2 x 2 + y 2 
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5.2 A final remark 



Thus, for elliptic equations of the second order with real-analytic coefficients 
in M 2 , there is no point to point reflection with respect to a real-analytic 
curve T unless T is a line or the following constrain a 2 + b 2 — 4c = for the 
coefficients of the equation holds. 

As it follows from |2QJ for elliptic equations in IR 2 , there is no point to 
finite set reflection as well. 

Point to compact set reflection is always possible. This set is a curve 
having one of its endpoints on a reflecting curve. The other endpoint is 
located at the reflected point itself. 
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